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Propagation of subcycle pulses in a two-level medium: Area-theorem breakdown and 

pulse shape 
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We solve the problem of ultrashort pulse propagation in a two-level medium beyond the rotating- 
wave (RWA) and slowly-varying-envelope approximations. The method of solution is based on the 
Maxwell-Bloch equations represented in the form that allows one to switch between RWA and 
general (non-RWA) cases in the framework of a single numerical algorithm. Using this method, the 
effect of a subcycle pulse (containing less than a single period of field oscillations) on the two-level 
medium was analyzed. It is shown that for such short pulses, the clear breakdown of the area 
theorem occurs for the pulses of large enough area. Moreover, deviations from the area theorem 
appear to be strongly dependent on the pulse shape that cannot be observed for longer few-cycle 
pulses. 

PACS numbers: 42.65.Re, 42.65. Sf, 42.50.Md, 42.65.Pc 



I. INTRODUCTION 

As the light pulses produced with modern laser tech- 
niques become shorter, so that their duration becomes 
comparable with the period of optical field oscillation, 
the necessity of an adequate theoretical description of 
the propagation of such ultrashort pulses in different 
systems tends to be more and more obvious. One of 
the basic models of an optical medium is a two-level 
medium, the fundamental model of a resonantly absorb- 
ing medium. The case of pulses containing only a few 
cycles of field oscillations implies that the analysis of the 
pulse dynamics in the two-level medium should be carried 
out beyond the popular and standard approximations - 
the slowly-varying-envelope approximation (SVEA) and 
the rotating- wave approximation (RWA). The study of 
few-cycle pulse propagation under such generalized con- 
ditions (beyond SVEA and RWA) has been under way 
since the mid 1990s and has given a number of impor- 
tant results. For example, the main effects previously 
known, such as self-induced transparency (SIT), 2w soli- 
ton formation, and 4n pulse splitting, were reported to be 
valid for the few-cycle pulse though some additional fea- 
tures in the dynamics of the two-level medium were found 
as well PHI]- Among other results one can mention the 
spectral transformations due to the intrapulse four- wave 
mixing Q, the local-field effects on the few-cycle pulse 
propagation d, [f|, the generation of a single-cycle soli- 
ton in a subwavelength structure consisting of the two- 
level medium Q , and the creation of the quasisolitons in 
a waveguide-like resonantly absorbing nanostructure Q ■ 
The most recent achievements include the effects of the 
chirp @-[l0| and the so-called counter-rotating terms in 
the Bloch equations on femtosecond pulse propaga- 
tion. 



In this paper we study the validity of the area theorem 
for even shorter (subcycle) pulses. Previously Hughes 
[l2| discovered the breakdown of the area theorem for 
the few-cycle pulses of large area (27rn with n > 3), 
while for lower areas this theorem is still able to predict 
the profile change of the pulse accompanied by its split- 
ting [l^ . Later Tarasishin et al. Q reported that the 
half- and quarter-cycle 2ir pulses leave some small resid- 
ual excitation inside the two-level medium. Here we show 
that these deviations from the area theorem are much 
more noticeable at larger areas of the incident subcycle 
pulses and we trace appearance of the theorem break- 
down with a shortening of the pulse. Moreover, these 
deviations appear to be strongly dependent on the pulse 
shape. The effect of pulse form on the excitation prob- 
ability of the two-level system is known for nonresonant 
excitation (see, for example, the work by Conover [l4| 
and references therein). However, we consider strictly 
resonant pulses, their shape being important only for the 
number of cycles less than unity. 

The structure of the paper is as follows. In Section 
lHl we establish the Maxwell-Bloch equations to be nu- 
merically solved and give the main parameters of cal- 
culations. In Section IIII1 comparing our results with 
the results known from the literature, we prove that the 
method based on the equations stated in the previous 
section can be applied to simulate ultrashort pulse prop- 
agation beyond the RWA and SVEA. Section [TV] is de- 
voted to the study of subcycle pulse interaction with the 
two-level medium, namely to the issues of the area theo- 
rem breakdown and the influence of pulse shape. Finally, 
in Section [V] we give a short conclusion. 



II. MAIN EQUATIONS AND PARAMETERS 



* dvnovitsky@tut.by 



Light propagation in the two-level medium beyond the 
RWA and SVEA is given by the Maxwell-Bloch equations 
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as follows [|[lH: 
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where E is the electric held of a light wave, P12 the off- 
resonant density matrix element (atomic polarization), 
w = PTi—P\\ the inversion (population difference), luq the 
frequency of atomic resonance, p the dipole moment of 
the quantum transition, 71 and 72 the rates of relaxation 
of population and polarization, respectively, c the speed 
of light, and ft the Planck constant. Here the macroscopic 
polarization of the two-level medium is P = 2/iCRepi2 
with C as the concentration (density) of two-level atoms. 
The symbols Re and Im stand for taking of real and imag- 
inary parts, respectively. 

Our aim is to rewrite Eqs. ([T]) to ([3]) in such a manner 
that they would allow direct comparison of the calcula- 
tions conducted with and without the RWA in the frame- 
work of a single numerical algorithm. To reach this aim, 
we represent the electric field and atomic polarization as 
E = {Acxp[i(ujt — kz)] + c.c.}/2 and P12 = pexp[i(u>t — 
kz)], respectively, but the complex amplitudes A and p 
are not assumed to be slowly varying. Here ui is the cen- 
tral frequency of radiation, k — u>/c is the wavenumber, 
and c.c. stands for complex conjugated term. Introduc- 
ing dimensionless arguments r = ut and £ = kz and the 
dimensionless field amplitude Q = {p. /Too) A (normalized 
Rabi frequency), we come to the set of equations 
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where 8 = Au/w = (u>q — uj)/lu is the frequency detuning, 
7i 2 = 71.2/^ are the normalized relaxation rates, and 
e = lol/oj = Att ji 2 'C l 'Zhui is the dimensionless parameter 
of interaction between light and matter (or normalized 
Lorentz frequency) . Finally, the auxiliary two- valued co- 
efficient s marks the situation considered: s — corre- 
sponds to the RWA (absence of "rapidly rotating" terms), 
while s = 1 is related to the general case. In this paper 
we numerically solve Eqs. ((H) to ©, so that we have the 
possibility of switching between the general (non-RWA) 



and RWA cases by simply choosing the appropriate value 
of the single parameter. The numerical approach is es- 
sentially the same as in our previous publications [TH-tHl 
where the relatively long pulses were studied in the limit 
of the RWA (but not the SVEA in the wave equation). 
Therefore, we do not discuss the details of the method 
and refer the reader to those works. 

We adopt the following parameters of the medium and 
light throughout the paper: the relaxation rates 71 = 1 
and 72 = 10 ns _1 are large enough so that we are in 
the regime of coherent light-matter interaction; the de- 
tuning (5 = (exact resonance); the light wavelength 
A = 2irc/ui = 0.83 pm; and the strength of light-matter 
coupling u)l = 10 11 s _1 which is much less than the ra- 
diation frequency. These material parameters mean that 
one needs to take relatively long thicknesses of the two- 
level medium (L > > A) to observe the transformations of 
the few- or subcycle pulse envelope. Therefore we do not 
discuss here the effects of pulse profile changing consid- 
ered previously for the relatively stronger coupling con- 
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Moreover, in 



the strong coupling limit, when the peak Rabi frequency 
flow is comparable to or less than the Lorentz frequency, 
one needs to take into account the so-called local field ef- 
fects [l7|. In our research we can neglect them, since the 
opposite inequality takes place (flow » ljl). The esti- 
mation shows that a single-cycle 27r soliton in a collection 
of two-level atoms with dipole moments /i~1D should 
have the peak amplitude of about 0.4 GV/cm. The re- 
quired concentration in this case is C ~ 2 x 10 19 cm -3 . 
Obviously, the results of the calculations can be rescaled 
for another set of parameters without loss of generality. 

In this paper we consider the pulses of two different 
shapes: the hyperbolic secant = fl p sech(t/t p ) and 
Gaussian 17 = Q p exp(— t 2 j2t 2 v ). The duration of the 
pulse t p is defined through the number of cycles N as 
t p = NT/f, where T = A/c is the period of electric 
field oscillations, and / is the coefficient which depends 
on the pulse form and describes its full width at half 
maximum (FWHM). For the hyperbolic secant pulse we 
have / = 2 arccoshy / 2, while for the Gaussian shape 
/ = 2\/ln 2. The peak (normalized) Rabi frequency 
flp is measured in the units of fig corresponding to the 
area of 2ir, so that for the hyperbolic secant pulse one 
should take £Iq = \/irct p and for the Gaussian pulse 



Oo = A/v27rcip. The calculational region includes the 
two-level medium of thickness L surrounded by the vac- 
uum regions of length 0.64 /im. The medium is supposed 
to be initially in the ground state (w = —1). 



III. TESTING THE CALCULATION 
APPROACH 

First of all, we need to ascertain that the method based 
on solving of Eqs. ^ to ^ correctly describes propaga- 
tion of the light pulses in the two-level medium. There 
are two such tests that are to be considered further. 
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FIG. 1. (a) The profile of the Gaussian 50-cycle 3tt pulse 
transmitted through the two-level medium of length L — 
100A. Calculations were carried out for the general case 
(s = 1). (b) The difference between the intensity profiles 
obtained for s — 1 and s = 0. (c) The difference between the 
intensity profiles calculated by two different RWA schemes: 
the scheme of Eqs. Q to © at s = and that of Ref. [l]| . 
Intensities are normalized by 7o = 0%. 
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FIG. 2. (a) Dynamics of the inversion at the entrance of 
the two-level medium excited by the Gaussian two-cycle 2n 
pulse. Calculations were carried out both for the general case 
(s = 1) and for the RWA (s = 0). (b) The difference between 
the inversion profiles obtained for s — 1 and s = 0. (c) 
The difference between the inversion profiles calculated by 
two different RWA schemes: the scheme of Eqs. Q to ([6]) at 
s = and that of Ref. [TJ]. 



(i) The limit of long pulses. In this situation one ex- 
pects that the calculations at s — 1 and s = give 
approximately the same result. To prove these expec- 
tations, we launch the Gaussian 37r pulse of 50 cycles 
(the duration t p « 83 fs) into the two-level medium of 
thickness L = 100A = 83 /im. We also carry out the 
calculations according to the numerical scheme of Ref. 
fl6j to check the consistency with the RWA case realized 
in our previous works fla - [l9{ . Figure QTa) shows the in- 
tensity profile of such a long pulse transmitted through 
the layer of thickness L = 100A; calculations were per- 
formed by Eqs. Q to ^ at s = 1. One of the main fea- 
tures of such a long pulse dynamics is seen: the pulse is 
compressed while forming the constant- form soliton [l8[ . 
Simulations of the RWA scheme show the agreement with 
the profile of Fig. [Ha). The plots in Figs. QJb) andQJc) 
demonstrate the accuracy of this agreement: the differ- 
ence between the intensity profiles obtained for s = 1 
and s = is as low as several hundredth of io (the unit 
of intensity equal to 0§ ) , while the discrepancy between 



the calculations by two RWA schemes does not exceed 
10 _4 Jo, respectively. This proves the correctness of our 
approach in the long pulse limit. 

(ii) Behavior of inversion in the case of a few-cycle 
pulse. Since the shape of the few-cycle pulse varies too 
slowly as it propagates, our second test deals with the 
time variation of the inversion w at the entrance of the 
two-level medium. To excite the medium, the two-cycle 
2tt Gaussian pulse is used (its duration is about 3.32 fs). 
The results of calculations are demonstrated in Fig. [3] 
The difference between s = 1 and s — cases is clearly 
seen and reaches values as large as 0.15. The inversion 
profile calculated without the RWA shows also the feature 
characteristic for the few-cycle pulse propagation - the 
step-like flattenings corresponding to the extremes of the 
time derivative of the electric field [l[. Finally, Fig. [5Jc) 
is the evidence of the precise correspondence between the 
calculation at s = and the RWA calculations according 
to the previously used approach. Thus, we conclude that 
our numerical method based on Eqs. ^ to ^ allows one 
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FIG. 3. (Color online) (a) Dynamics of the inversion at the 
entrance of the two-level medium excited by the Gaussian and 
hyperbolic secant 3n pulses of different numbers of cycles N. 
(b) The corresponding dependence of the final state of inver- 
sion (FSI) Wf on the number of cycles. FSI was determined 
at the time point t = 50t p . The inset shows the field profiles 
of the half-cycle 3n pulses of the Gaussian and sech shapes. 

to reproduce the main peculiarities of long- and short- 
pulse propagation discovered previously and can be used 
to study the subcycle pulses beyond the RWA. 

IV. RESULTS ON SUBCYCLE PULSES 

Let us apply the method described above to simulate 
propagation of subcycle pulses in the two-level medium. 
The main parameter to be traced is the final state of in- 
version (FSI) denoted here by Wf. It is the steady value 
of inversion in which the medium appears after passage 
of the incident pulse. The area theorem implies that 
FSI depends on the pulse area: if the starting state of 
inversion is w s = —1 (two-level system in the ground 
state), then for the pulse area mr one will have Wt = 1 
or —1 at odd and even values of integer n. Figure [SJa) 




FIG. 4. (Color online) The dependence of the final state of 
inversion (FSI) Wf on the area of the half-cycle pulses of the 
Gaussian and hyperbolic secant shapes. The results for the 
RWA case are given for comparison. FSI was determined at 
the time point t = 50t p . 

shows the dynamics of inversion at the entrance of the 
two-level medium under the influence of the 3n subcycle 
pulse, i.e. for the number of cycles N < 1. It is seen that 
for the single-cycle pulse (N — 1) the FSI approximately 
(though not exactly) corresponds to the value predicted 
by the area theorem (wf = 1). As the duration of the 
pulse decreases (N — 0.75 and 0.5), the breakdown of 
the area theorem becomes apparent: such a short pulse 
simply cannot guarantee the full cycle of inversion dy- 
namics, so that w finally appears somewhere in between 
— 1 and 1. The concrete value of Wf strongly depends on 
the pulse shape. We considered two different variants - 
secant hyperbolic and Gaussian pulses - and see that the 
difference between the FSI in these two variants grows 
as the pulses become shorter. For comparison, we also 
calculated the RWA curves for the half-cycle pulses [see 
the lower right panel in Fig. [Ha)]. In this case the area 
theorem is strictly valid for the pulses of any shape. In 
other words, the importance of the rapidly rotating terms 
in the Bloch equations increases for ultrashort pulses and 
results not only in the breakdown of the area theorem, 
but also in the strong dependence of the medium dynam- 
ics on the pulse form. This is seen in the dependence of 
Wf on the number of cycles N shown in Fig. EJb): for 
the Gaussian subcycle pulses, Wf deviates faster from 
the area theorem than in the case of sech pulses. This 
result seems to be rather surprising, since the field pro- 
files of sech and Gaussian half-cycle pulses seem to be not 
so much different (see the inset) to lead to such a large 
difference in Wf (Awj w 0.35). Thus, in the range of sub- 
cycle pulses, even a small distinction in the pulse shape 
may result in a strong change of the medium dynamics. 
This should be taken into account when performing ex- 
periments with subcycle pulses and gives an additional 
parameter to control the state of the medium. 

The next issue is the dependence of the above described 
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breakdown of the area theorem on the pulse area. Let 
us calculate the FSI for different areas of the incident 
half-cycle pulse (N = 0.5). The results for the Gaussian 
and hyperbolic secant forms are presented in Fig. 2] as 
well as the RWA data which are the same for both pulse 
profiles. It is seen that at small areas all three curves 
approximately coincide and only small deviations from 
the area theorem occur as was reported previously [3|. 
But at the areas above 2n the curves rapidly diverge, so 
that the large-area half-cycle pulses carry the medium to 
the state with the level population which depends on the 
shape of the pulse. Moreover, the positions of minima 
and maxima of inversion also shift from the usual values. 
This is especially clearly seen for the second minimum: 
while it is the area of 4tt in the RWA case, it is signif- 
icantly lower in the general case - about 3.9-7T for the 
hyperbolic secant pulse and even 3.7-7T for the Gaussian 
one. It is also worth noting that the curve for the sech 
shape is closer to the RWA dependence in the vicinity of 
37r areas, but the situation turns out to be opposite in 
the vicinity of 5ir where the curve for the Gaussian pulse 
approaches the values consistent with the area theorem. 
The need of large areas (> 2ir) to observe the strong de- 
viations from the area theorem can be explained by the 
fact that the medium affected by the subcycle pulse does 
not have enough time to develop the full cycle of dynam- 



ics when it includes more than a single excitation and 
deexcitation. 

V. CONCLUSION 

In summary, we have implemented and tested the nu- 
merical approach which allows one to solve the Maxwell- 
Bloch equations beyond the RWA and SVEA. This 
method was used to study the effects of the interaction of 
the subcycle pulses with the two-level medium. For such 
pulses containing less than one period of optical field, we 
have found that (i) the breakdown of the area theorem 
becomes apparent at areas larger than 2tt, and (ii) the 
result of light-matter interaction strongly depends on the 
pulse shape. Our research was focused on the study of 
the medium dynamics so that we have not considered the 
effects of subcycle pulse profile change such as the forma- 
tion of an asymmetric wave form |3j. Another interest- 
ing issue worth to be studying further is the possibility or 
impossibility of subcycle self-induced-transparency (SIT) 
solitons. On the one hand, there are studies that show 
the existence of few-cycle SIT solitons 20] . On the other 
hand, the attempts to find a solitonic regime for subcy- 
cle pulses were not successful so far [3]. The approach 
reported in this paper can be used for a detailed study 
of this question. 
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